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Abstract 



It is widely believed that some classes of high-energy transients may be 
powered by the rotational energy of a rapidly spinning black hole. The energy 
extraction mechanism commonly discussed involves macroscopic magnetic fields 
that are produced by currents flowing in a disk or torus surrounding the black 
^T) . hole. The discovery of relativistic jets in radio loud AGNs nearly half a century 

O^l ' ago was the main motivation for the analysis presented in the seminal paper 

by Blandford and Znajeck, where a global model for the magnetosphere of a 
Kerr black hole in the force-free limit was first constructed and employed to 
demonstrate that, under reasonable astrophysical conditions, energy can be 
extracted quite efficiently in the form of a Poynting flux. The recent discoveries 
of relativistic motions in Galactic microquasars, and the indications that GRBs 



eject ultra-relativistic, collimated outflows, have lent support to the hypothesis 



that this mechanism may be universal. 



In spite of large efforts that have led to impressive progress in our under- 
standing of the physics of magnetized flows in Kerr geometry, several impor- 
tant issues remain unresolved: the nature of the load in the Blandford-Znajeck 
, model, the causality and stability of a force free magnetosphere, global current 

in ' closure, and the role of boundary conditions on the channels through which 

the extracted energy is released. This chapter provides an account of recent 
work that addresses some of those open questions. After a brief introduction, 
an overview of MHD in Kerr geometry will be given. The basic properties of 
a static black hole magnetosphere will be discussed, with particular empha- 
sis on the conditions under which efficient energy extraction may occur. The 
trans-field equation that determines the geometry of magnetic flux surfaces, 
and the role of its critical surfaces will be summarized. The force-free limit 
will be considered, with an attempt to clarify several issues concerning the long 
standing causality problem. This will be followed by a discussion on global 
magnetosphcric structures and related emission channels. A particular model, 
developed recently for the interaction of a uniformly magnetized torus with a 
rapidly spinning black hole under extreme conditions, those anticipated in mag- 
netars and GRBs, will be considered in some greater detail. The magnetic field 
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configuration invoked in this model is vastly different from the standard one 
adopted in previous work, and leads to interesting astrophysical consequences 
that will be described in detail. In particular, the model predicts that the major 
fraction of the black hole spin energy should be released in the form of gravi- 
tational waves, with estimated detection rate of up to a few per year with the 
LIGO and VIRGO detectors. Non-stationary effects will be considered in the 
last part of this chapter. The requirements for local adjustments of a magnetic 
flux tube to temporal changes will be carefully examined, and it will be argued 
that changes in the angular velocity of magnetic lines induced by small distur- 
bances, involve effects beyond the lowest order geometric optical approximation. 
The global evolution of a force- free magnetosphere will be explored next, with 
an attempt to elucidate the critical role played by the frame-dragging dynamo. 
Recent analytic results will be presented and contrasted with full GRMHD 
simulations. 



1 Introduction 



It is widely believed that some classes of high-energy transients, notably AGNs, 
microquasars, and gamma-ray bursts (GRBs), may be powered by the rotational 
energy of a rapidly spinning black hole. The energy extraction mechanism com- 
monly discussed involves macroscopic magnetic fields that are produced by currents 
flowing in a disk or torus surrounding the black hole. A version of this idea was 
proposed originally by Ruffini k Wilson (1975). Blandford k Znajek (1977; hence- 
forth BZ) were the first to construct a global model for the magnetosphere of a Kerr 
black hole in the force- free limit, and to demonstrate that energy can be extracted 
under certain conditions in the form of a Poynting flux. They applied their model to 
radio loud AGNs, and proposed that extragalactic jets are the consequence of this 
energy extraction mechanism (see review by Begelman et al. 1984). The dissipation 
of the Poynting flux is not addressed within the framework of the BZ model; it is 
merely assumed to occur in some non-force-free region where the magnetic field is 
sufficiently weak. The detection of prodigious gamma-ray emission from blazars by 
EGRET on CGRO (Thompson et al. 1993) indicates that a considerable fraction of 
the bulk energy is dissipated already on rather small scales. Conversion of magnetic 
energy to kinetic energy is a feature that appears to be common to many classes 
of relativistic astrophysical systems. The mechanism responsible for this process is 
poorly understood at present, although some ideas have been discussed (Romanova 
k Lovelace 1992, 1997; Levinson 1996, 1998; Thomson 1997; Levinson k van Put- 
ten 1997; Kirk k Lyubarsky 2001; Lyubarsky 2003). An important diagnostic of 
the BZ model is the composition of the collimated outflows, which is expected to 
be dominated by an electron-positron plasma. Unfortunately, the content of extra- 
galactic as well as Galactic jets is yet an open issue (e.g., Mannheim 1993; Dermer 
k Shlickheiser 1993; Sikora et al. 1994; Blandford k Levinson 1995; Celotti 1997). 
Some constraints can be imposed from the observations (Sikora et al. 1997; Celotti 
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1997), but the analysis is inconclusive. This issue may be eventually resolved by the 
upcoming km 3 neutrino detectors. Various aspects of the BZ model have been subse- 
quently examined in more detail by different authors (e.g., Bekenstein & Oron 1978; 
Macdonald & Thorne 1982; Phinney 1983; Thorne et al. 1986; Punsly & Coroniti, 
1990; Ghosh & Abramovitz 1997), but many issues remained to be resolved. 

The BZ model has been applied in recent years also to stellar systems, in partic- 
ular GRBs (e.g., Levinson & Eichler 1993; Meszaros & Rees 1997; Van Putten 2000, 
2001; Lee, et al. 2000; Brown, et al. 2000), and microquasars (Levinson & Blandford 
1996). The formation of systems consisting of a rapidly rotating, stellar mass black 
hole surrounded by a magnetized torus, is thought to be the outcome of catastrophic 
events such as black hole- neutron star and neutron star-neutron star coalescence 
(e.g., Eichler et al. 1989; Paczyhski 1991) or core collapse of massive stars (Woosley 
1993). The putative existence of such black hole-torus remnants is by itself a strong 
motivation to consider the physics and observational consequences of the interac- 
tion of the black hole and the torus. It is quite likely that upcoming gamma-ray 
(e.g., MAGIC, GLAST), km 3 neutrino (e.g., Ice-Cube, NEMO), and gravitational 
wave (LIGO, VIRGO) detectors will detect new signatures in currently identified 
classes of high-energy transients, or may even discover new classes of astrophysi- 
cal objects, that may be the product of the underlying black hole-torus systems. 
A particular example is emission of gravitational waves by a torus in suspended 
accretion state (Van Putten 2001, 2002; Van Putten & Levinson 2002, 2003). A 
key feature in this model that distinguishes it from the original BZ scenario, is a 
strong coupling between the black hole and the torus which is mediated by angular 
momentum transfer along magnetic field lines connecting the horizon and the inner 
face of the torus. The emission of the gravitational waves by the torus, as well as 
intense baryon rich winds is a consequence of this coupling. This example suggests 
that remnants of collapsed, massive stars may be prodigious sources of gravitational 
waves. The proposed association of these objects with presently known classes of 
GRBs (Van Putten 2001, Van Putten & Levinson 2001), is motivated in this model 
both phenomenologically and from theoretical considerations that imply that a small 
fraction of the black hole rotational energy should be expelled along the axis in the 
form of ultra-relativistic, baryon poor jets. However, the properties of the emission 
from those baryon poor jets are likely to depend also on environmental conditions 
and additional microphysics, and it could well be that the association with iden- 
tified classes of GRBs applies only to a sub-class of these potential gravitational 
wave sources. Moreover, the prompt GRB emission is presumably strongly beamed, 
while the gravitational wave emission is roughly isotropic, so that most of those pu- 
tative gravitational wave sources are not expected to be coincide with an observed 
GRB anyhow. On the other hand, delayed, isotropic radio emission is anticipated 
from the interaction of the torus winds with the ambient medium, and/or from the 
GRB afterglow. Systematic searches for radio transients may eventually reveal the- 
ses remnants (Levinson et al. 2002). A search for associations of radio transients 
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with LIGO signals may provide valuable information and may aid detecting these 
sources. 

The total amount of energy released in the BZ process depends mainly on the 
mass, M, and angular momentum, a, of the black hole, and the corresponding rate 
depends on the strength of the poloidal magnetic field threading the horizon, B. To 
be more precise, the Blandford-Znajek power can be expressed as (see section |2~5|) : 

PB Z = e[^) 2 B\ 2 H c, (1) 

where rn = 2GM/c 2 is the Schwartchild radius, and e is an efficiency factor that 
depends on the geometry of the magnetic field, as will be shown below. Typically, e 
lies in the range 10~ 2 — 10 _1 . For a supermassive black hole of mass M ~ 10 9 solar 
masses, as expected in the powerful blazars, magnetic field of the order of 10 4 G is 
required to explain the energetics of the associated gamma-ray jets. The magnetic 
field may be supported by a surrounding disk and extend close to the event horizon 
(e.g., Bisnovatyi-Kogan & Ruzmaikin 1976). In the case of GRBs, that harbor 
stellar mass black holes, super-critical magnetic fields, of order 10 15 G are needed 
in order to account for the characteristic luminosities measured. Such fields are 
much stronger than those inferred in typical radio pulsars. Which process can give 
rise to such amplification of the magnetic field is yet an open issue, although some 
models have been proposed (e.g., Thomson & Duncan 1993). Recent observational 
efforts provide indications of supercritical magnetic fields in magnetars (Duncan & 
Thompson 1992; Kouveliotou et al. 2003), suggesting that such high magnetic fields 
may also be present in other extreme systems. 

2 IDEAL MHD IN KERR GEOMETRY 

2.1 Basic equations and integrals of motion 

We express the Kerr metric (Kerr 1963) in Boyer-Lindquest coordinates with the 
following notation: 

n 2 

ds 2 = -a 2 dt 2 + u 2 (dcf) + (idtf + ^-dr 2 + p 2 d6 2 , (2) 

where a = p\fK/Yi is the lapse function, Co 2 = {Y? / p 2 ) sin 2 9, and —(3 = 2aMr/T, 2 
is the angular velocity of a ZAMO with respect to a distant observer , with A = 
r 2 + a 2 — 2Mr, p 2 = r 2 + a 2 cos 2 9, and T? = (r 2 + a 2 ) 2 — a 2 A sin 2 9. The parameters 
M and a are the mass and angular momentum per unit mass of the black hole. 

We denote by n, p, p, h = (p + p)/n, respectively, the proper particle density, 
pressure, energy density, and specific enthalpy of the MHD flow. The stress-energy 
tensor then takes the form: 

T Q/3 = hnu a vP + pg aP + ^(F a(7 F^ - \g a(3 F 2 ), (3) 

An 4 
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where u a is the four- velocity, and F^ v = d^A^ — d v A^ is the electromagnetic tensor. 
The dynamics of the MHD system is then governed by the following set of equations: 
Maxwell's equations, 

Ff a a = J=(V=p^) )a = 4vr/, (4) 

+ Fp~t,<* + -^7",/3 = ^' (5) 



the continuity equation 



K) ;a = -^=(V=gnu a ) >a = 0, (6) 



and the energy and momentum equations 

T^ v = -±=(^Tn,v + KfiT * = °- (7) 

Here denotes the associated Cristofel symbol 

In regions where the MHD flow is ideal (i.e., infinite conductivity), Ohms law 
yields the additional constraint, 

u a F afS = 0. (8) 

The stationary axisymmetric flow considered here is characterized by two Killing 
vectors: £^ = dt and = 9<j>- By contracting these Killing vectors with the stress- 
energy tensor we can construct the energy and angular momentum currents, which 
in Boyer-Lindquest coordinates read: 

F = = -T' t = -hnU r U t - ^F r9 F te , (9) 

£ e = T 6 e = -T t 9 = -hnU e U t + ^-F r0 F tr , (10) 

C = 7>V = Tl = hnU r U s + ^F^g, (11) 

C e = 7>V = Tl = hnU 9 ^ - ^F^F+r. (12) 

These currents are conserved, viz., 

£% = E* = 0. (13) 

The relation 

= FtfFfr + Fr+Ftf = -(F^gdr + F ri pdo)A ( f > = D^As, (14) 

where the operator = Fgsdr + F^dg denotes derivative along magnetic flux 
surfaces, implies that As is conserved on magnetic flux surfaces and is, therefore, a 
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viable stream function. We can use the dual electromagnetic tensor, denoted here 
by , to express the invariant E • B in the form, 

E-B = -(V=ff/4)^F^ = £)^ t . (15) 

It can be readily shown that the ideal MHD condition, given by eq. (JEJ), implies that 
E • B = 0, meaning that the electric potential is conserved along magnetic surfaces. 
We note that in cases where a gauge At = —QA^ can be found, where 0(r, 0) is some 
function of the coordinates, the change in the electric potential along flux surfaces 
is given by D^A t = —A^D^Q, where eq. (fTI"]) has been used. Thus, in regions 
where the MHD flow is ideal f2 is conserved along magnetic surfaces. As discussed 
below the function f2 is approximately the angular velocity of the magnetic flux tube 
at small angles. On the other hand, in regions where the ideal MHD condition is 
violated, e.g., the sparking gap (see below), the change in angular velocity across 
this region is related to the electric potential drop through: AAt = — A^Atl (note 
that in steady-state A^ is conserved even in the non-ideal case). The relevancy of 
this result to the double transonic flow is discussed below. 

Equations (JSJ), ©, ©> and (|§|)- ((T3|) admit four additional quantities that are 
conserved along magnetic flux surfaces (e.g., Bekenstein &: Oron 1978): The particle 
flux per unit magnetic flux, 



ata V-gnu ^-gnu 

vm = — ^ — = — p — ; (16) 

the angular velocity of magnetic field lines, 

n F (V) = - jJL j Fro + v*; (17) 
y — gnu 1 

and the total energy and angular momentum per particle carried by the MHD flow, 



E(9) = -hu t ~ ^-n F F re , (18) 

47T7/ 

L(9) = hu, - £1f«. (19) 

Here = u^/u 1 , and v r = u r /u t are the corresponding components of the 3- 
velocity. The fifth integral of motion is the entropy s = s(^f). Equations (|16|) - (|19|) 
can be solved for Uf, u§ and F r6 . The solution can then be used to express Fg, 



a' = g u ut + g^Uff,, and u^ = g^ut + g^u^ as: 



* — y "t T y 1 «0, aim w — y at T y 1 1 

2 r ,~,2 



47TT? a 2 L-^^l F + (3){E + ^L) 
e ~ sm9 a 2 -Cj 2 {VL f + (i) 2 - M 2 ' 1 ' 

. t _ 1 a 2 (E-n F L)-M 2 (E + f3L) 

a 2 a 2 -Co 2 {VL F + (5) 2 -M 2 ' 1 j 

_ a 2 Co 2 n F (E - n F L) + u 2 !3M 2 {E + (3L) - a 2 M 2 L 
~ a 2 ti 2 [a 2 -ti 2 (Q F + f3) 2 - M 2 } ' [ ' 
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where M 2 = {^ithrj 2 ) jn. The meaning of the latter becomes clear when written in 
terms of ZAMO 4- velocities. By employing eq. ((TB|) one obtains M 2 = a 2 u 2 /u\, 
where u p = (u r u r + -u% e ) 1/2 and u A = (B^/Anhn) 1 / 2 , with B 2 = (Esinfl)" 2 ^ + 
AF 2 r ), are, respectively, the poloidal velocity and Alfven velocity, as measured by a 
ZAMO (see e.g., Macdonald & Thorne 1982, Beskin 1997). Thus, M is up to a factor 
a the Alfven Mach number. Finally, using the normalization condition, u^u^ = — 1 
and the definition of the Mach number M, we obtain the Bernoulli equation 

M 4 B 2 

U " + l= IforW* 2 + 1 = {aUt? " ^ W + U0)2 ' (23) 
Note that since the enthalpy h is a function of the thermodynamic variables n 
and s only, the density and enthalpy can be expressed in terms of s and M, viz., 
n = n(s,M), h = h(s,M). Thus, equations l|20j ) -(|2$j ) determine essentially all the 
flow characteristics in terms of the five integrals of motion once the poloidal magnetic 
field is known. 



2.2 Asymptotic behavior of the MHD flow 

Consider first the behavior of the solution near the horizon. There a — > 0, and for 
physical solutions for which the Mach number is finite on the horizon, M(th) 7^ 0, 
eqs. (|21j) and (|22|) yield — > Sl#u* as a — > 0, where VLh = —j3n is the angular 
velocity of the black hole. Consequently, the rotation of the plasma on the horizon is 
synchronous with the black hole, as one might expect. Substituting the latter result 
into Bernoulli equation (|23|) . we obtain u p — > au l on the horizon. The poloidal 
velocity is radial on the horizon, implying u 2 — ► g rr u r u r , hence 

v r = u T /u t ^a/^=-A/(r 2 + a 2 ). (24) 

This shows that the plasma moves along geodesies of a freely falling observer as it 
approaches the horizon, meaning that near the horizon the dynamics is governed by 
gravity alone. Substituting and v r obtained above into eqs. (fTo^l and (|17j). yields 

^ = - r -^(Q F -Q H ). (25) 

Eq. (|25j) gives the frozen-in condition derived originally by Znajek (1977) and used 
by BZ as a boundary condition in their force-free analysis. 

Next, consider the behavior of a MHD outflow far from the black hole. The 
outflow parameters are given, to a good approximation, by eqs (|20 |) -([23 |) with a 2 = 1, 
P = 0. At the Alfven surface the Mach number is M\ = 1 — £b 2 A Q 2 F: and the 
requirement that Fg u l and remain finite there yields the well known relation 
L/E = uj^^lp (e.g., Weber & Davis 1967; Camenzind 1986). By employing eqs. 
(|TT)|) and (fTT)) we obtain, 

F re M 2 (R 2 nE - L) {Cj 2 - Co 2 A )M 2 £l F 



F^g R 2 v r [(l - M 2 )E - m] v r Co 2 {M 2 - M\) 



(26) 
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Well above the Alfven point, where Co » Co a and M >> Ma the last equation gives 

^ - — • (27) 

The asymptotic poloidal current follows from the r component of eq. (|4"|): 

<9 / A sin 



= 47T v /=?j r . (28) 



50 V p 2 

Integrating the latter equation we obtain the net electric current within a flux tube, 

A sin 6* 
V 

Using eqs Q25|l and (|27() we find that the net electric current flowing on the horizon 
is given by 



2n^fde = -Z^-Fro. (29) 



S sin 

and the net electric current at infinity is given by, 



lH = ^^(^H-n F )F^ e , (30) 



Ioo = Y^sm.eF <i>6 , (31) 

2.3 Conditions for energy extraction by a magnetized inflow 

Detailed analysis of the inflow structure and the requirements for extraction of the 
hole rotational energy is given in Takahashi et al. (1990) and Hirotani et al. (1992). 
For completeness, we give in this section a brief derivation of the conditions under 
which energy extraction by a MHD inflow is possible. From equations (|9" ]) -(|12 )) . (|18|) 
and (j!9j) it is readily seen that the energy and angular momentum fluxes can be 
expressed as: E a = Enu a , and C a = Lnu a , with a = r,6. Here nu a is the particle 
flux carried by the MHD flow. Consequently, energy extraction by the MHD inflow 
(for which nu a is negative) then requires the specific energy to be negative, viz., 
E(ijj) < 0. This shows that the BZ mechanism is indeed a manifestation of the 
Penrose process, as pointed out by Takahashi et al. (1990) 

Now, at the Alfven surface the denominator of eqs ()20|) - (|22j) vanishes, and we 
obtain 

Ml = a 2 A -Cj A {n F + p A ) 2 > 0, (32) 

where the subscript A refers to quantities on the Alfven surface. This in turn implies 
that the angular velocity of magnetic field lines must lie in the range 

- (3 A - a A /CoA < < ~Pa + oa/Coa- (33) 
In order for to be finite on the Alfven surface, the relation 

L _ Cb 2 A (Q F + p A ) 



E a A -Q A P A (n F + /3 A ) 



(34) 
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must hold (see eq. |21j). At the injection surface of the inflow we take M 2 ~ 0. 
Equation (|2*3*)) then yields 

E-Q F L = h inj [a 2 inj - u 2 nj (Q F + /3 inj ) 2 ] 1/2 > 0. (35) 

Solving eqs (|3*1|) and (|3*5|) we obtain 



_ a 2 A -U 2 A P A (^F + M h r 2 ~2 / 0e , , „ ^1/2 

- — 2 ^TTF, . n n 2 inj [inj inj\ F ' Ptnj) \ , (3bj 

a A ~ ^AV^ ~r PA) 

T — & a (SIf + Pa) , r 2 _ ~2 i o \2ll/2 /o 7 \ 

a A~ u A." l F + PA) 



Evidently, the specific energy will be negative provided the condition oP a —Cj a (5a{S^f+ 
(3a) < is satisfied on the Alfven surface. The latter condition combined with 
eq. Q33|) implies, in turn, that the rotational energy of the hole can be extracted 
only if (i) the Alfven surface of the inflow is located inside the ergosphere, and (ii) 
< VL F < Vt H . 



2.4 The trans-field equation 

Equation (fT6|) can be used to express the poloidal components of the 4-velocity in 
terms of the magnetic field components. Substituting the latter into eq. (jSJ, the 
poloidal components of the stress-energy tensor can be written as: 

M 2 1 1 

T ab = i6?r2 2 (d a if>)(M) - P9 ab + ^{F aa F? + -g^F 2 ), (38) 

where ifi = ItxA^ is the magnetic flux, and the indices a, b run through r and 8 only. 
The poloidal magnetic and electric fields in eq. (|38|) are given in terms of the stream 
function ijj as: F^ a = — (l/27r)9 a V> F ta = —QpF^a, and the toroidal magnetic field 
F r g by eq. (|20[). Upon substituting eq. (|3*5|) into eq. (J7J), one obtains a second 
order PDE for the stream function ip, involving the five integrals of motion, the 
Mach number, the proper density n and the temperature. The resultant trans-field 
equation is given explicitly in eq. (39) of Beskin (1997) 1 . Combined with Bernoulli 
equation (|23|) and an appropriate equation of state, that allow one to express, though 
implicitly, the Mach number and the thermodynamic parameters in terms of the 
conserved quantities, this trans- field equation can be solved in principle to yield the 
stream function, i/j = ip(r,9), for a given set of boundary conditions. The functional 
form of the invariants , L(ifj), E{tp), r](ip) and s(ip) that appear explicitly in 
the equation must be specified first through appropriate boundary conditions. 

The trans-field equation has, in general, several singular surfaces. Those singular 
surfaces can be identified from differentiation of eq. I)23j) . whereby an expression for 

1 The trans-field equation is written in Beskin 1997, who uses the 3+1 formalism, in terms of the 
operator V = e r \fKp~ 1 d r + egp~ 1 dg 
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the gradient of the Mach number can be obtained (Takahashi 1990; Beskin 1997). Of 
most interest are the Alfven surface 2 , on which the Mach number satisfies relation 
(|32|) and the fast and slow magnetosonic surfaces on which the poloidal 4-velocity 
satisfies 

Up = K± ^K 2 - \c 2 u\ [1 - (fl F + /3) 2 w 2 a- 2 ], (39) 

with the + (— ) sign applies to the fast (slow) magnetosonic surface. Here 2K = 
[1 - (fi F + (3) 2 u 2 a- 2 ]u 2 A + (B 2 /47r/m) + c 2 , where = (^/p 2 )F rd is the ZAMO 
toridal magnetic field, and c s is the sound 4-velocity, given by c 2 = a 2 /(l — a 2 ) with 
a 2 = (l/h)(dp/dn) s . As the flow accelerates, it crosses first the slow surface, then 
the Alfven surface, next the light cylinder, and finally the fast surface. As shown in 
e.g., Beskin (1997), with the exception of the force-free case the trans-field equation 
changes from elliptic to hyperbolic on the fast critical surface, meaning that the flow 
beyond that surface (between the fast critical surface and the horizon in the inflow 
section) is not in causal contact with the region between the fast critical surface 
and the injection point. In the force- free limit discussed next, the Alfven surface 
coincides with the ligh cylinder, and the fast critical surface with the horizon. In 
this case the trans-field equation remains elliptic down to the horizon. Now, the 
requirement that the MHD flow pass smoothly through the critical surface imposes 
certain relations between the conserved quantities of the MHD flow. Since the trans- 
field equation is second order the number of boundary conditions required is 2 + 
number of conserved quantities - number of singular surfaces. 

There is ample literature on the properties of the Grad-Shafranov equation and 
its solutions in flat spacetime, mainly in the context of pulsar physics (see e.g., 
Mestel 1999, and references therein). Solutions of the trans-field equation in curved 
spacetime were obtained in various limits. Examples are discussed in, e.g., (Bland- 
ford & Znajek 1977; Fendt 1997; Beskin 1997; Ghosh 2000; Beskin & Kuznetsova 
2000; Lee, at al. 2001; Uzdensky 2004a,b) 

2.5 The force-free limit 

In situations where the MHD flow is so highly magnetized that the inertia of the 
plasma can be ignored, the system approaches the force-free limit. Formally this 
limit can be obtained by setting h — > in the above equations 3 . This in turn 
implies M 2 — > (meaning essentially that the Alfven velocity, as measured by a 
ZAMO, approaches the speed of light). The Alfven surface then approaches the 
light cylinder, as can be readily seen from eq. ()32j) . In the limit h = eq. (J7J) 

2 It is noteworthy that although the Grad-Shafranov equation has a singularity on the Alfven 
surface, the regularity of equations 12111 - 122(1 is automatically satisfied. It merely defines the location 
of this surface. 

3 We note that in reality h > mc 2 , where m is the rest mass of a particle in the flow. The 
limit M 2 « 1 is due to the small particle density. Thus, the limit of negligible inertia is better 
characterized by the dimensionless Mach number. 
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reduces to 



F^f = 0, 



(40) 



where eq. Q has been used. The toroidal component of the last equation and the 
condition F^ t = give F$ r f + F^f = 0, meaning that the poloidal current must 
flow along magnetic flux surfaces. Consequently, the net electric current within a 
given flux surface, given by eq. (|29jl. is conserved, as can be readily verified by using 
the r and 9 components of eq. and the latter condition. 

The total energy and angular momentum are given by Etot = TjE and Ltot = ijL. 
By employing eqs. (|T8|) and (|T9|) we obtain E to t = ^F-^tot (to be more precise, eq. 
|35| shows that this result is correct, in general, to order hi n jrj). Taking the limit 
M 2 — >■ and E — » VlpL in eq. (|20j) and using eq. (j29|) we also find / = lixLtot- 
Thus, the force-free system has only two integrals of motion, Q,p and I. 

Combining eqs © and ()25j) . one can obtain the energy flux on the horizon: 



which is a well known result, derived originally by BZ. Maximum extraction effi- 
ciency occurs for Vip = The net power transferred outwards from the horizon 
is P = 2ir J ^/—g£ r d0, where the integration is over the horizon surface. By express- 
ing F^g in terms of the ZAMO poloidal field, B r = —F^fEsmO, in eq. (|41|). and 
scaling the metric components, we arrive at eq. Q, with the dimensionless factor e 
given explicitly as an integral involving the scaled components. 

2.6 The Double transonic flow 

In many applications of the BZ mechanism to astrophysical systems it is envisioned 
that the energy extracted near the horizon of the hole is ultimately transferred to 
large radii in the form of a relativistic jet. Since the poloidal velocity of the out- 
flowing plasma is in the opposite direction to that of the negative energy inflow 
inside the ergosphere, the invariants r](ip), E{^i) and L(ip) must have an opposite 
sign in the two regions (the energy and angular momentum fluxes can be continues 
nonetheless), implying that there must exist a region, that we shall term the in- 
jection region, where the ideal MHD condition is violated, and where fresh plasma 
is produced. This region must be located somewhere between the inner and outer 
Alfven surfaces (Takahashi et al. 1990). We point out that this discontinuity in 
the flow does not appear implicitly in the force-free case, where formally r] = 0. As 
shown above, the force-free system has only two invariants, / and both of which 
can be assumed continuous in the entire region between the inner and outer fast 
magnetosonic surfaces. The implicit assumption being made is that there exists a 
source of electric charges, presumably located on the neutral surface defined below. 

For illustration let us consider the properties of the double transonic flow near 
the rotation axis of the black hole (figure 1). We suppose that the MHD inflow 



s r = n F (n H -n F )-(F (t>e ) 2 



(41) 
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Figure 1: Schematic illustration of a double-transonic structure. Transonics inflow and out- 
flow are expelled from the injection region where plasma is continuously generated. Magnetic 
flux surfaces are assumed to be continuous across the injection region. The angular veloci- 
ties of magnetic field lines in the inflow and outflow sections are determined by the potential 
drop AV along magnetic surfaces in the injection region, and the cross-field electric current 
J s = I H - loo (see text). 
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and outflow expelled from the injection region become ideal outside the injection 
region. The entire system is then characterized by 10 invariants, 5 for each flow. We 
henceforth refer to quantities in the inflow (outflow) by a subscript + (-). Consider 
first the charge distribution near the axis. The t component of eq. (j3J) reduces in 
the small angle approximation to, 



d_ 

de 



or 



47TA 



-93 



(42) 



For reasonable magnetic field configurations eq. (|42|) can be integrated to yield the 
lowest order Goldriech-Julian charge density: 



(Sl F + P)B r 
2vra 2 



(43) 



where B r is the ZAMO radial magnetic field. As seen, the charge density changes 
sign across a neutral sheet on which the angular velocity of magnetic lines, as mea- 
sured by a ZAMO is zero, viz., Q,p = —(3. Now, to lowest order eqs. ()3Uj) and 1)31(1 
give 



(fin - n+)ip 



2tt 



2tt 



(44) 
(45) 



Assuming that the stream function is continuous across the injection region then 
implies that an amount 



AI = I H - I 



2tt 



(46) 



of cross-field electric current must be supplied into the flux tube. In situations where 
inertia is negligible, so that the inflow and outflow can be considered force-free, the 
current AI must be supplied by the injection region. Under the assumption that 
the stream function is continuous a gauge can be found such that At — —^/2tt to 
lowest order. The potential drop along magnetic flux tubes in the injection region 
is then given by (see eq ^2] and the text below), 



AV 



- o_)^/27r. 



(47) 



Once a solution ij)(r, 0) of the Grad-Shafranov equation is found, eqs. (|46|) and (I47|) 
can be solved to yield the angular velocity of magnetic field lines on the horizon 
near the axis in terms of the magnetic flux, cross field current, and potential drop 
in the injection region: 



2vr(A/ - AU) 

4> 



(48) 
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In general, AI and AV will be determined by the conditions in the injection region. 
Consequently, the angular velocity of magnetic lines on the horizon and, hence, 
the efficiency at which energy is extracted from the hole are determined by the 
microphysics of the injection region. The values of AV and AI, the toroidal electric 
current in the interface separating the inflow and outflow, and the continuity of 
the stream function across this interface provide 4 boundary conditions for the 
Grad-Shafranov equation governing the structure of the double transonic flow. The 
values of the injected densities, rii n j±, enthalpies , hi n j±, and poloidal velocities, 
u p ±\inj, provide the six additional boundary conditions required. Thus, it seems 
that the entire structure is determined essentially by the physical conditions in the 
injection region and the regularity conditions at the critical surfaces alone. No 
additional boundary conditions are required on the horizon. This remains true in 
the limit of very small (albeit finite) inertia. In particular, when AI << Iqj and 
AV « V max , as expected in situations whereby the plasma source is associated 
with a quasi-steady sparking gap, inside which the field aligned electric field is 
almost completely screened out, we obtain £1 + ~ 0^/2 near the axis, confirming 
earlier results (Blandford & Znajek 1977, Phinney 1983). Detailed analysis of a 
perturbed split monopole solution for a double transonic flow is presented in Beskin 
& Kuznetsova (2000), who reached similar conclusions. 

2.7 Inflow from a uniformly magnetized, massive torus 

As a second example we consider a situation in which the magnetic field lines extend- 
ing from the horizon are anchored to a massive torus. This is shown schematically 
in figure 2. We suppose that the mass flux along field lines connecting the torus 
and the black hole is small, and approximate the inflow as force-free. This may 
be justified in cases where angular momentum transfer from the hole to the torus 
results in a slightly super-Keplerian rotation of the inner parts of the torus, thereby 
strongly suppressing accretion (van Putten & Ostriker 2001). 

For an infinitely conducting torus, the angular velocity of magnetic field lines 
that originate from the torus equals the the angular velocity of the torus, viz., 

= 0^. The torque exerted on the inner face of the torus by the black hole is 
obtained by integrating the angular momentum flux C r given by eq. (|11|) over the 
section of the horizon which is threaded by magnetic field lines that are anchored 
to the torus: 



where eq. (|25|) has bee used. Here Oh is the angle of the last field line that connect 
the torus and the horizon. In terms of the net poloidal magnetic flux associated 
with the open field lines in the torus, ipo, we may write r = (Tin — ^t) fn^o / ^ • 
Formally fn is defined through eq. (|49j) . It roughly represents the fraction of surface 
area of the horizon that is threaded by magnetic field lines emerging from the torus. 




(49) 
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Figure 2: Inflow from a massive, infinitely conducting torus: In this example there is a 
direct link between the torus and the black hole. The angular velocity of magnetic field lines 
that penetrate the horizon equals fir- Energy and angular momentum are transferred from 
the black hole to the torus when fin > £^t, and in the opposite direction when fin < fir- 
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It is seen that when £Ih > Qt angular momentum is transferred from the black 
hole to the torus, tending to spin up the inner face, whereas in the slowly rotating 
case (Qh < ^t) the black hole receives angular momentum from the torus. The 
associated power is given simply by L = Qtt. 

2.8 Is there a causality problem? 

The question whether a force-free magnetosphere can exist and whether it is stable 
has been the subject of a long debate (Punsly & Coroniti 1990, Beskin & Kuznetsova 
2000; Komissarov 2001,2003,2004b; Blandford 2001, 2002; Punsly 2003, 2004 van 
Putten & Levinson 2003; Levinson 2004). In particular it has been argued that a 
force-free black hole magnetosphere is not a causal structure and, therefore, physi- 
cally excluded. The claim made by Punsly and collaborators is that the fact that the 
angular velocity, £If, of a force-free flux tube extending from the horizon is not a free 
parameter, but is determined by matching boundary conditions on the horizon and 
at infinity (as in the analysis of BZ and Phinney 1983) violates the principle of MHD 
causality. The reason is that the inflowing magnetic wind must pass through the 
inner light cylinder before reaching the horizon and, therefore, cannot communicate 
with the plasma source region (see section H~Tj) . Their main conclusion is that the 
use of the Znajek frozen-in condition on the horizon to determine £lp is unphysical, 
and that Qf must be determined by the dissipative process that leads to ejection 
of plasma on magnetic field lines between the inner and outer Alfven points. This 
view seems to follow from the interpretation of the membrane paradigm, that the 
horizon can be viewed as a conducting surface (see also Komissarov 2004b). It has 
already been emphasized in the preceding section that the Grad-Shafranov equation 
is in general hyperbolic near the horizon and requires no boundary conditions there, 
and that the Znajeck frozen-in condition is essentially a regularity condition on the 
fast critical surface. Moreover, it has been argued that the injection region alone 
determines the entire solution of the double transonic flow. The question whether a 
force- free magnetosphere is stable and how it evolves is nonetheless relevant. Bland- 
ford (2001, 2002) suggested that any changes are communicated by means of a fast 
magnetosonic mode, that can propagate across magnetic field lines at the speed of 
light and carry information about the toroidal magnetic field and poloidal current 
to the plasma source even beyond the light cylinder. Punsly (2003) argued that the 
fast mode can only carry displacement currents, and therefore concluded that fast 
characteristics cannot transmit sufficient information to affect the global structure. 
Levinson (2004) in turn claimed that to second order the fast mode can affect the 
poloidal current (see section ETTl below for further details). The present view of this 
author, as explained below, is that although this is an interesting issue, it is of little 
relevancy to the causality problem. 

To elucidate how a force-free magnetosphere responds to changes in spacetime, 
consider the following Gedunken experiment. Suppose that a double transonic flow 
has been ejected along the rotation axis of a Kerr black hole, and assume further that 
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the flow is force-free everywhere except for the plasma source region where a small 
deviation from force- freeness must exist. Imagine now that the angular momentum 
of the hole changes abruptly at time to from ao to ao + Sa. This will induce a change 
—5(3{r) in the angular velocity of a ZAMO at any radius r. The question then is 
how this variation in space-time would affect the flow. From eq. (|4S[) it is evident 
that a change 



in the local electric charge density near the axis is required in order for the MHD 
flow to remain ideal, where 8£If is the resultant change in the angular velocity of 
the magnetic flux tube. This is particularly true in the region between the inner 
Alfven and fast critical surfaces. Whether the charge density adjusts by means 
of second order effects of the fast mode, as suggested by Levinson (2004; but cf. 
Punsly 2004), or otherwise is an interesting question by its own. In our view it 
must adjust. The crucial point is that any deviation of j* from the Goldriech-Julian 
value will induce a field-aligned electric field (i.e., E • B ^ 0), in violation of the 
force-free assumption, that would tend to be screened out instantaneously by the 
surrounding plasma. Such a process involves most likely generation of longitudinal 
modes that cannot be analyzed within the framework of force-free wave analysis 
anyhow. Whether a quasi steady-state can be maintained in such a case is unclear 
to this author. Recent time-dependent analysis of sparking gaps demonstrates that a 
response of this kind might be oscillatory in nature (Levinson et al., in preparation). 
The analysis outlined in section l4~2l demonstrates that the magnetosphere remains 
force-free during the course of evolution provided that sufficient plasma is present 
on field linse. 

Now, if conditions in the plasma source are such that it always tends to be sus- 
tained in an approximately force-free state, in the sense that the cross field current 
AI and the potential drop across the plasma source region AV are negligibly small, 
then it follows from eqs. (jHOj) and IjHlj) that as long as Qf / + <5^#)/2, the 
charge density in the gap will change with time sine V • j 7^ 0. This will feedback 
on the potential drop AV, and will lead to an evolution of the poloidal current and, 
hence, Qf, until the value (Qh + 8Qh)/2 is exceeded 4 . 



The global structure of the magnetosphere determines the geometry of magnetic 
flux surfaces and the location of regions where substantial departure from ideal 
MHD takes place. The existence of such regions is essential in order to allow for (i) 
cross-field currents that are necessary for global current closure, and (ii) dissipation 
of magnetic energy to account for the high entropy inferred in most astrophysical 

4 This assumes again that the system can relax to a steady-state. Alternatively, it may undergoes 
large amplitude oscillations. 




(50) 
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Figure 3: Schematic representation of an open-field magnetosphere 

applications. Moreover, the global structure determines through which channels the 
BH rotational energy extracted is emitted. 

In most astrophysical applications it is envisioned that the magnetic field thread- 
ing the horizon is supported by currents flowing in a disk or torus surrounding the 
black hole. Two inherently different disk magnetizations are discussed in the liter- 
ature. The first one consists of what we term an open-field magnetosphere. In the 
second one, termed closed-field geometry, a large portion of magnetic field lines that 
thread the horizon are anchored to the torus (Nitta et al. 1991; van Putten 1999, 
2001; Uzdensky 2004b). This introduces a strong coupling between the black hole 
and the torus. In what follows we discuss those two magnetospheric configurations 
in some greater detail. 



3.1 Open-Field Geometry 

This configuration has been studied extensively (e.g., Blandford 1976; Lovelace 1976; 
Phinney 1983; Macdonald 1984; Beskin & Parev 1993), and was the one employed 
originally by BZ in their model. It is shown schematically in fig 3. In this config- 
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uration there is no direct link between the horizon and the disk. In particular, all 
field lines that penetrate the horizon extend to infinity. The energy extracted from 
the hole is then transported outward to the load along those field lines in the form 
of a Poynting flux. The nature of the load, where magnetic energy is presumably 
being converted to kinetic energy, is not well understood at present. The current 
closure path is also not well understood. It is conceivable that a current sheet may 
form on the equatorial plan, thereby providing a return path. This idea seems to be 
supported by recent numerical simulations (Komissarov 2004b). If so, then a torque 
may be exerted on the inner parts of the disk by the cross-field currents. This may 
give rise to enhanced dissipation, which has claimed to be inferred recently from 
X-ray observations of MCG-6-30-15 (Wilms, et al. 2001). However, rapid heating 
of the inner parts of the disk are expected also in the closed-field magnetosphere, 
though for a different reason as explained below. 

3.2 Closed-Field Geometry 

The basic properties of such a magnetosphere can be illustrated by matching a 
vacuum solution to a uniformly magnetized torus. Such a solution is topologically 
equivalent to a configuration produced by two counter-oriented current rings in the 
equatorial plane, which in flat space-time can be calculated analytically. The resul- 
tant magnetic flux surfaces of this vacuo solution are exhibited in fig. 4. A third 
current loop associated with the black hole equilibrium magnetic moment in its 
lowest energy stat (van Putten 2001) has been added in the figure. It is oriented 
antiparallel to the surrounding torus magnetosphere, facilitating an essentially uni- 
form and maximal horizon flux at arbitrary spin-rates. It only affects the solution 
very near the horizon. As seen, there are two separated regions of closed field lines. 
At large radii (compared with the radius of the outer current ring) the field quickly 
approaches a dipole solution. In the inner region the field lines intersect the hori- 
zon, giving rise to a strong coupling between the black hole and the inner face of 
the torus. 

Now, the magnetic field inside the torus cannot be purely poloidal, because 
purely poloidal fields are unstable and tend to decay completely in a few Alfven 
timescales (Markey k Tayler, 1973; Flowers & Ruderman, 1975; Eichler 1982). 
However, by conservation of helicity, a twisted magnetic field does not decay to zero, 
at least in the limit of ideal magnetohydrodynamics. Instead, it will evolve into a 
new, stable configuration. Recent 3D MHD simulations (Braithwaite & Spruit, 2004) 
show that magnetic fields inside stars tend to develope a belt of twisted field lines 
that stabilize a dipolar field in the magnetosphere above the stellar surface. This 
configuration appears to be stable over the resistive timescale, which is typically 
much longer than the canonical dynamical timescales (e.g., rotation periods and 
accoustic timescales). By topological equivalence in poloidal cross-section, it is 
therefore conceivable that the magnetic field inside the torus is twisted as well, 
supporting an overall torus magnetosphere which, from the outside, is consistent 
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Figure 4: The poloidal topology of magnetic flux-surfaces in vacuum produced by two 
counter-oriented current rings is shown. The curent rings are located on the equatorial 
plan, one on the side of the torus facing the black hole and the other one on the opposite 
side, as indicated. The dashed line is the separatrix between the flux-surfaces supported by 
the inner and the outer faces of the torus. Reprinted from van Putten (2004). 
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Figure 5: Schematic illustration of the poloidal topology of the magnctosphere of a torus 
surrounding a rapidly rotating black hole. The direct link between the inner torus magne- 
tosphere and the horizon results in energy and angular momentum transfer from the black 
hole to the torus. This input is catalyzed by the torus into gravitational radiation, winds, 
thermal and MeV neutrino emissions. The associated Maxwell stresses are mediated by 
poloidal currents, which close over a current sheet (marked as vacuum gap) in an annulus of 
vanishing magnetic field (as indicated). A baryon-poor inner flux-tube serves as an artery 
for a minor fraction of black hole-spin energy. The inner and outer flux tubes are separated 
by a charge and current sheet in the outflow section. The dashed lines indicate the inner 
and outer light cylinders. The lifetime of the system is set by the lifetime of rapid spin of 
the black hole. Reprinted from van Putten & Levinson ApJ 584, 937 (2003). 
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with a uniformly magnetized surface of the torus. 

Like in pulsars (Goldriech & Julian 1969; Michel 1982), the vacuum magneto- 
sphere of the rotating torus is unstable. By vacuum break-down, the flux-surfaces 
will evolve with electric charges to a largely force- free state. As a result, a mag- 
netosphere develops which consists of conductive flux-surfaces and magnetic winds. 
This is shown schematically in fig 5. In the limit of infinite torus conductivity, the 
flux-surfaces in the outer /inner torus magnetosphere assume rigid corotation with 
the outer/inner face of the torus. A transonic inflow is expelled from the inner face 
along the open field lines to the horizon. For a rapidly rotating black hole we antici- 
pate < and so energy and angular momentum are transferred from the hole 
into the inner face of the torus (see section 2.7), tending to spin it up. Likewise, 
a transonic outflow is ejected from the outer face to infinity, resulting in a loss of 
energy and angular momentum of the outer face. A quasi steady-state is quickly 
reached, whereby a flow of angular momentum mediated by shear forces between 
the differentially rotating torus layers is established (van Putten & Ostriker 2001). 
The energy deposited in the torus is emitted predominantly in the form of gravita- 
tional waves, owing to large deformations of the torus, as well as MeV neutrinos and 
baryon rich winds to infinity, resulting from the rapid heating of the torus by the 
friction between its layers (van Putten & Levinson 2003). This is one of the major 
differences between the open-and-closed- field geometries. 

When the magnetic field strength inside the torus exceeds a few time 10 15 G 
it becomes dynamically important. It is then anticipated that magnetic stresses 
that builds up inside the torus will lead to nonlinear, dynamical deformations of the 
torus. The power spectrum of the induced mass moments is likely to be dominated 
by the several lowest multipoles. A naive estimate of the field strength above which 
the torus becomes unstable to deformations may be obtained by equating the torque 
acting on a perturb current ring, owing to the mutual magnetic interaction between 
the two current rings, with the gravitational torque exerted by the central black hole 
(van Putten & Levinson 2003). This yields a critical field strength of about 10 16 G. 
The presence of toroidal field components inside the torus may somewhat alter this 
estimate. Interestingly, this critical magnetic field corresponds to a spin down time of 
the order of tens of seconds for a stellar mass black hole, consistent with durations 
of long GRBs. Full 3D GRMHD simulations are ultimately needed to study the 
dynamics of the torus (including the back-reaction of the black hole) in this regime. 
Using a toy model, Bromberg et al. (2005) have recently calculated the evolution of 
the mass moments of the deformed torus. Their simulations reveal the existence of 
a nonlinear, oscillatory phase that sets in when the magnetic field approaches the 
critical value. The duration of this phase is found to be very long compared with the 
orbital period, for a certain range of parameters. The deformations of the torus that 
result from this magnetic self-interaction should give rise to a burst of gravitational 
wave emission with durations of several to several tens of seconds (see Bromberg et 
al. 2005 for analysis of the gravitational wave spectrum). For reasonable parameters 
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we anticipate such sources to be detectable by LIGO and VIRGO out to a distance 
of 100 Mpc and even larger (van Putten et al. 2004) 

The powerful wind driven by the pressure gradients in the surface layers of the 
hot torus passes through the outer Alfven point. This results in the opening of some 
magnetic field lines in the outer layers of the inner torus magnetosphere. Because the 
torus is rotating and magnetized, the ejection of the wind is partially anisotropic. 
Specifically, mass flux is generally suppressed along magnetic field lines that are 
inclined toward the rotation axis, and enhanced along field lines that are strongly 
inclined away from the axis (Blandford &; Payne 1982; Romanova et al. 1997), 
owing to centrifugal forces. The details of the outflow depend on the heating and 
cooling rate of the corona, and on its structure. We speculate that large pressure 
gradients in the corona would tend to push matter along some of the magnetic field 
lines originally connected to the horizon, and that a combination of buoyancy and 
centrifugal forces may subsequently give rise to a twist of these field lines, some 
of which may ultimately fold and open to form an open magnetic flux tube near 
the rotation axis that extends from the horizon to infinity. The resulting structure 
consists of two coaxial flux tubes with opposite magnetic orientation, as shown in 
fig 5. The inner and outer flux tubes are separated by a cylindrical current and 
charge sheet that accounts for the jump in the electric and magnetic fields across 
the interface. The lower section of the inner/outer flux-tube which connects to 
the horizon has, instead, a parallel orientation between the poloidal magnetic field. 
In the perfect MHD limit, the properties of the interface are described by jump 
conditions as follow from Maxwell's equations. The surface charge density is given 

by 

4vro- e = [QrsmOB r ] = rs\n0tt T B T+ - r sin6»(0^/2)5 r _, (51) 

where B r+ (-Br-) denotes the radial magnetic field near the interface in the outer 
(inner) flux tube, and the poloidal and toroidal surface currents by, 

4vrJ r =[B^ = 
4vrJ^ =-[B r ]. 

The poloidal current (|52|) results beyond the Alfven point, where the wind transports 
angular momentum outwards to infinity. The outflow in the inner tube is expected 
to be relativistic {v r = 1). If the baryon rich outflow from the torus also becomes 
relativistic, then the latter equation implies that J r = cr e c, namely the poloidal 
current in the boundary layer is solely due to the outflowing surface charges. This 
current sheet is a potential site for reconnection of magnetic field lines, which would 
convert magnetic energy in the inner flux tube into kinetic energy. 

4 NON-STATIONARY EFFECTS 

The discussion in the preceding sections focused on static magnetospheres. The 
question of whether these structures are globally stable and how they evolve can 
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only be addressed using time dependent analysis. This is the ultimate goal of general 
relativistic MHD (GRMHD) simulations, that only very recently became feasible. 
Simple analytic analysis is nonetheless important, both to give us some insight into 
this complicated problem and to provide test cases for the simulations whenever 
possible. 

There are two crucial issues that we think need to be clarified. The first one is 
the role of linear waves in the response of a magnetized flux tube to local changes. 
The second and most important one, is the role frame dragging plays in the global 
evolution of a force- free magnetosphere. We discuss these issues in tern in what 
follows. 

4.1 Small amplitude waves 

The system of GRMHD equations admits in general 4 types of modes - slow, Alfven, 
fast and entropy, that can most easily be identified from a linear perturbation anal- 
ysis in a homogeneous background. There is ample literature on the properties of 
those waves (for recent accounts see e.g., Anile 1989; Uchida 1997a,b; Komissarov 
1999, 2002; Punsly 2001; van Putten 2004). The linear modes also define the criti- 
cal surfaces of a static transonic flow as discussed in section l2~4l above. In essence, 
each such surface acts as a one-way membrane for waves of the corresponding type. 
Whenever the flow is locally perturbed, only information associated with subcritical 
modes can be transmitted backwards to the source. In particular, no information 
can be transmitted from the regions beyond the inner and outer fast critical sur- 
faces, meaning that the horizon of a Kerr black hole is not in causal contact with 
the injection region. The question then arises, how information about the state of 
the black hole is transmitted through a magnetized flux tube. 

In the force-free limit only two physical modes remain, the Alfven and fast modes 
(e.g., Uchida 1997b; Komissarov 2002). As explained in section T2.5( in this limit 
the Alfven surface coincides with the light cylinder, whereas the fast critical surface 
approaches the horizon. This fact has led Blandford (2002) to propose that changes 
near the horizon can be communicated by means of the fast mode, a proposal that 
has been questioned by Punsly (2001; 2003). Below, we re-examine this problem. 
The most detailed investigation of force-free waves in Kerr spacetime is presented 
in Uchida (1997a, b). This author derived a set of PDEs for the Lagrangian dis- 
placement of linear perturbations. Using a WKB approximation he then solved it 
to lowest order. His results confirm that to lowest order the fast wave is electromag- 
netic, in the sense that it propagates along null geodesies and is purely transverse. 
The Alfven mode has, in general, a longitudinal component. It also propagates at 
the speed of light, but only along poloidal magnetic field lines. It is these proper- 
ties of the force-free waves (see Punsly 2003) that are at the base of the causality 
dispute discussed in section |2~%1 As we shall now argue, the response of a flux tube 
to changes in spacetime involves higher order effects, which are not accounted for 
in the lowest order geometric optical approximation. To illustrate this, consider a 
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static flux tube extending from the horizon of a Kerr black hole. Using eq. I|43|) 
the GJ charge density near the axis, as measured by a ZAMO, can be expressed in 
terms of the 8 component of the ZAMO electric field as, 

PGJ = (53) 

where u is the cylindrical radius, as defined below eq. ©. Suppose now that the 
system is transformed from this initial state to a slightly different state, e.g., due to 
a change in the angular momentum of the black hole by a small amount 5a « a. 
The change in the angular velocity of the flux tube, as measured by a ZAMO, will 
be accompanied by a change 5Eq in the electric field, and a corresponding change 
in the GJ charge density: 

fe = si- < 54 > 

Consequently, the perturbed GJ charge density is associated with changes on scales 
comparable to the dimension of the flux tube. This means that the adjustment of 
the magnetosphere to global changes cannot be analyzed within the framework of 
the geometric optical approximation, as attempted by Punsly. It is not even clear 
whether for such long wavelength perturbations decomposition into MHD modes 
is possible. To be more precise, the perturbed charge density associated with any 
short-wavelength disturbance of some static solution is given by 5p e = zk • <5E/47r, 
where k is the corresponding wave vector. Thus, any changes in the charge density 
required for adjustments of the angular velocity of an evolving magnetic flux tube 
would appear only to order (/ccD) -1 , which is neglected in the lowest order geometric 
optical approximation. The fact that the fast mode is purely transverse to lowest 
order does not by itself imply that a force-free magnetosphere cannot respond to 
changes beyond the inner light cylinder. It simply means that account of higher 
order terms in the linear perturbation analysis is mandatory for exploring such 
effects. Such an attempt is presented in (Levinson 2004), who confirms that to 
lowest order the fast mode is indeed electromagnetic. However, he finds (but cf. 
Punsly 2004) that to second order the fast mode have a longitudinal component, 
and that the perturbed electric charge density beyond the light cylinder approaches 
$PGJ, as given in eq. (|54|). 

4.2 The frame dragging dynamo 

To gain some insight into the role of frame dragging in the evolution of a force- 
free flux tube, let us first express Faraday law in terms of the ZAMO fields Bs, 
E r , and Eg, and the potential As, which will be used as our free variables. In 
Boyer-Lindquist coordinates eq © gives, 



B<f,,t + -^-{paE e ) )T - -^{paE r ) t e = -^((3 tr As,e - 0,6As, r ). (55) 
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The term on the R.H.S of the last equation, which is absent in flat spacetime, repre- 
sents the effect of frame dragging. As seen it only couples to the poloidal magnetic 
field, and can be interpreted as a driver or a dynamo term. This term cannot be 
gauged away owing to the fact that frame dragging is differential. Physically, this 
can be attributed to the fact that the angular velocity of a ZAMO depends on ra- 
dius, and so magnetic surfaces appear to be in differential rotation in the ZAMO 
frame, thereby giving rise to a potential drop along magnetic surfaces that tends to 
be screened out when sufficient plasma is present. Let us explore further the role of 
this driving term in the evolution of a flux tube. 

Consider the evolution of a force-free magnetosphere of a slowly rotating black 
hole. To derive the equations for the evolving electromagnetic field, we linearize 
Maxwell equations Q and using the hole angular momentum a as the small 
parameter. We suppose that initially the magnetosphere is non-rotating, and that 
the magnetic field can be described by the vacuum Wald solution (Wald 1974). To 
second order in a/M the initial solution reads: 

A$ = (B/2)r 2 sm 2 9, (56) 

= E r = Eg = 0. To this order the force-free condition (jID*|l reduces to 



Aj r = 
AEg 



-r cot 9jg, 
r cot 8E r 



(57) 
(58) 



where j a denotes the components of the ZAMO poloidal current, which are related 
to the Boyer-Linduist current through: j r = (p/y/A)j r , jg = pj e . The r and 9 
components of eq. (jU yield, 



E r ,t + 



n 



psm 

a 

E$ t H — 
P 



(sin BB^fi = 4iraj r , 
= -Airaje, 



(59) 
(60) 



From equations (|55|) - (|6Ujl . we obtain a differential equation for the electric field Eg 
(see Levinson 2004 for further details): 



1 + A tan 2 6/r 2 )Egj t + 



AB d 



,t),r 



r 3 cos 6 



{sm.eB^ t ) fi = 0, 



with 



B, 



<j>,t 



AEg 



+ 



A, QMBa 

(tan 0Eg)fi ^ — sin a cos i 



(61) 



(62) 



The last term on the R.H.S of eq. 1)621) accounts for the differential frame dragging 
of the initial magnetic field, that is, P^A^^g — f3pA§^ r = — (6MBa/r 2 ) sin 6 cos 9 to 
second order in a/M. Now, since the magnetosphere is initially non-rotating, no 
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electric current is flowing in the system, viz., j r (t = 0) = jg{t = 0) = 0. To examine 
how the poloidal current is generated, we take the time derivative of eq. Q59|) . and 
employ equations and (j62j) to obtain near the rotation axis (that is, at small 
angles), 



v A 

j rt (t = 0) ~ -12BMa— r cos 2 0. (63) 

Consequently, the poloidal current is driven solely by the frame dragging dynamo; 
the only assumption being made is that there is sufficient plasma in space to al- 
low the condition E • B = to be satisfied everywhere. Note that the poloidal 
current is generated initially everywhere in space and not only on the horizon, in 
contrast to the case of a Faraday disk, as in Punsly's (2001) waveguide model (cf. 
Komissarov 2003). This clearly shows that the interpretation of the horizon as a 
unipolar inductor is inappropriate; it is the gravitomagnetic effect that is responsi- 
ble for the adjustment of the magnetosphere to changes. A similar conclusion was 
drawn earlier by Komissarov (2003, 2004b). Equation is rather complicated. 
Approximate, analytic solutions can be obtained near the horizon (Levinson 2004) 5 . 
For the case considered here we obtain in the region where a « 1, and for times 
t < 2Mln(2M 2 /A): 

E e = |^\/Asm20[cosh(i/2M) - 1], (64) 

B^ = -|^\/Asin20sinh(V2M), (65) 
313(1 

j r = - ^3 VA cos 2 9 sinh(t/2M). (66) 

As seen, the perturbed force-free field grows exponentially with an e-folding time 
2M. Note that after time t ~ 2Mln(2M 2 /A), at which the above solution is no 
longer valid, the toroidal magnetic field and the current evolve close to their steady- 
state values, viz., Bt ~ (Qh /tyF^g sin 6>, I = Bt/2. We conjecture that after this 
time the system will reach a steady state. This suggests that the magnetosphere 
inside the ergosphere evolves to a steady-state solution over a few dynamical times. 



4.3 Numerical Simulations 

There have been several attempts in recent years to perform numerical simulations 
of the Blandford-Znajek process. They all show the tendency of the magnetosphere 
to evolve towards a stable steady state. In the GRMHD simulations reported in 
Koide et al. (2002) and Koide (2003), the initial magnetic field configuration is 
described by the Wald solution. The plasma around the black hole has initially zero 
momentum, a uniform mass density, and Alfven velocity of 0.983 c (Koide 2003). 

5 There is a typo in eq. (24) of Levinson (2004). It should read: Bt, t = 
x2M 2 sin 9 (f - ff|| sin 26>) . 
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This initial condition is similar to the one invoked in section H~2l except that plasma 
inertia is not neglected and the black hole in the simulations is nearly maximally 
rotating. The simulations clearly show the generation of a toroidal magnetic field by 
the frame dragging dynamo inside the ergosphere, and the consequent decrease of the 
specific energy of the plasma until it becomes negative. The twist of the magnetic 
field lines seems to be propagating outward, carrying energy and angular momentum 
on account of the black hole rotational energy. The simulations reported in Koide 
et al. (2002) run for a time of about several r s /c after which the code crashes. This 
run is, unfortunately, not sufficiently long for the system to reach a steady-state. 
Nonetheless, it does show that the magnetosphere inside the ergosphere evolves 
over a few dynamical times, consistent with the analytic result derived in section 
14.21 Komissarov (2001) performed 2D, time dependent numerical simulations of 
the evolution of a force-free, axisymmetric monopole configuration around a Kerr 
black hole. He found that the solution quickly settles to a stable steady state, with 
the angular velocity of magnetic field lines approaching f2#/2 (meaning maximum 
extraction efficiency), thereby confirming the results of BZ. He later generalized the 
numerical model to include inertial terms (Komissarov 2004a), and demonstrated 
that the inertia of the plasma near the event horizon is dynamically unimportant, 
and that the force-free limit is a good approximation for magnetically dominated 
flows, at least for the monopole case. The simulations also show the development of 
a double transonic flow that extends beyond the outer fast critical surface. In a later 
paper (Komissarov 2004b), this author explored more realistic configurations. He 
also included a prescription for electric resistivity that enabled him to incorporate 
dissipative magnetospheric regions (current sheets) in the numerical model. He 
confirms his earlier conclusions, that a stable steady state is reached whereby energy 
is extracted electromagnetically. He finds, however, that in certain magnetospheric 
configurations, current sheets may form on the equatorial plane. Such regions may 
provide a path for the return current of the global circuit. Moreover, in certain cases 
collimation of the outflow is not seen up to a few tens of gravitational radii. 

In another set of numerical experiments (Hirose et al. 2004; McKinney & Gam- 
mie, 2004; De Villiers et al. 2005), the evolution of a weakly magnetized torus 
surrounding a Kerr black hole has been examined. It is generally found that the 
magnetic field in the torus is initially amplified via the MRI and becomes turbulent. 
A funnel region near the rotation axis of the hole is identified in those simulations, 
in which the magnetic field appears to be ordered and nearly force-free. This region 
is well described by the BZ model ( McKinney & Gammie, 2004), and it appears 
that energy is being extracted along those force-free flux tubes. However, the overall 
energy flux (integrated over the horizon) is dominated by the enthalpy of accreted 
plasma. This is not surprising in view of the initial and boundary conditions in- 
voked. Nonetheless, the presence of a magnetically dominated polar region along 
which a significant fraction of the accreted energy is channeled into a Poynting flux 
jet is interesting and of direct astrophysical relevance. 
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The conclusion to be drawn from the numerical experiments described above is 
that the magnetosphere of a Kerr black hole is a stable, causal structure, and that 
its evolution is driven by the frame dragging dynamo. 

5 CONCLUSION 

• The system of static, axisymmetric ideal MHD equations in Kerr geometry 
is characterized by five quantities conserved along magnetic flux surfaces; the 
specific energy and angular momentum, the angular velocity of magnetic lines, 
the ratio of particle and magnetic fluxes, and the entropy. The stream function 
that defines those flux surfaces obeys a second order, nonlinear PDE that 
involves the five invariants. 

• The specific energy of a given magnetic surface that extends down to the 
horizon becomes negative when two conditions are satisfied: (i) its angular 
velocity is larger than zero and smaller than that of the black hole, and (ii) 
the corresponding Alfven point is located inside the ergosphere. Energy and 
angular momentum can be transmitted from the horizon outwards along such a 
negative energy flux tube, on account of the black hole rotational energy. The 
rate at which energy is extracted from the hole depends on the strength of the 
magnetic field near the inner fast critical surface, and the angular velocity of 
the flux tube. For typical astrophysical parameters, the total power that can be 
extracted is sufficient to account for the luminosities exhibited by relativistic 
systems, e.g., AGNs, GRBs, microquasars, provided the efficiency is high. 

• Transonic solutions must pass through several singular surfaces. Of most in- 
terest are the slow magnetosonic, Alfven and fast magnetosonic surfaces, on 
which the bulk velocity of the flow equals the velocity of the corresponding 
linear mode. There are two sets of such critical surfaces, the inner one associ- 
ated with a transonic inflow into the hole, and an outer one associated with a 
transonic outflow to infinity. Each surface acts as a one-way membrane to the 
corresponding mode, and determines what information can be communicated 
backward to the injection region. The entire structure of the flow is deter- 
mined by appropriate boundary conditions and the regularity conditions on 
the critical surfaces. This is true also in the limit of zero inertia. 

• Double transonic structures that extend from the horizon to beyond the outer 
fast critical surface, as in the applications to astrophysical jets, must con- 
tain a region where the ideal MHD condition is violated. This region serves 
as a plasma source and must be located somewhere between the inner and 
outer Alfven surfaces. The global double-transonic flow structure is controlled 
entirely by the micro-physical conditions in the plasma source and the regu- 
larity conditions. In the force-free case the deviation from force-freeness in the 
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plasma source may be very small, and in any case is required only to ensure 
the continuity of electric current along magnetic field lines. In some configu- 
rations, as in the closed-field geometry discussed in section 13.21 a cylindrical 
current sheet may separate the inner, double transonic flow and the outer 
magnetosphere. 

• The extracted energy may be released through various channels, depending 
on the global structure of the magnetosphere. In open-field magnetospheres, 
where there is no direct link between the black hole and the surrounding disk, 
the extracted energy is ultimately channeled along the rotation axis in the 
form of a Poynting flux dominated outflow. In closed-field magnetospheres, 
a significant fraction of magnetic field lines that penetrate the horizon are 
anchored to the surrounding torus. In such configurations the major fraction 
of the extracted energy is transfered to the torus, and only about 0.1 % are 
channeled along the rotation axis. The energy deposited into the torus may 
result in strong gravitational wave emission and baryon rich winds. The ram 
pressure of the baryon rich wind that ensheath the inner jet can provide a 
means for collimating the inner jet. 

• The global evolution of a magnetosphere is governed by frame dragging. It is 
the fact that the ZAMO angular velocity varies with radius that causes the ap- 
pearance of field-aligned electric field in vacuum (or a starved magnetosphere 
more generally). In regions where sufficient plasma is present the parallel elec- 
tric field tends to be screened out. The adjustment of the electric charge to 
changes in spacetime is instantaneous essentially, and is dictated solely by the 
requirement that the condition E • B = is preserved in the course of evolu- 
tion. The term rhB • V/9 that appears in the homogeneous Maxwell equations, 
where B is the ZAMO poloidal magnetic field, -(5 is its angular velocity, and 
rh is a unit vector in the (j) direction, can be treated as a driving term that 
generates the toroidal magnetic field and poloidal electric current. 
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